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Abstract 

We perform a theoretical study on direct CP violation in in phase 

space around the intermediate states p^(770) and /o(500). The possible interference 
between the amplitudes corresponding to the two resonances is taken into account, and 
the relative strong phase of the two amplitudes is treated as a free parameter. Our 
analysis shows that by properly chosen the strong phase, both the CP violation strength 
and differential decay width accommodate to the experimental results. 


1 Introduction 

Charge-Parity {CP) violation has gain extensive attentions ever since its first discovery in 
systems in 1964 [1]. Within the Standard Model (SM), CP violation is originated 
from a complex phase in Cabibbo-Kobayashi-Maskawa (CKM) matrix, which describes the 
mixing of weak and mass eigenstates of quarks [2]. 

Although it is a small effect in general, CP violation can be relatively large in some decay 
channels of B and Bg mesons [3, 4, 5]. In fact, large CP violation has been confirmed in some 
two-body decay channel of B and Bg meson, such as B^ T^CP(+i)^^ [6], 

[7], B^ [8], B^ /0(1370)7r^[9], B^ p-K+ [10], and Bg tt+R- 

[11, 12]. In recent years, even larger CP violation which localized in three-body decay phase 
space of B meson was observed by LHCb collaboration in channels such as B^ 
and B^ TT [13, 14]. In view of its anisotropy property for small invariant mass 

of 7r+7r“ pair (no larger than the mass of p^(770)) in phase space, the large localized CP 
violation was first interpreted as a consequence of the interference of the decay amplitudes 
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corresponding to nearby resonances with different spins [15, 16]. Some other explanations 
such as the re-scattering effects of final states [17, 18] were also proposed thereafter. 

Since it is believed to be very small within SM, CP violation in the charm sector provides 
a good place for searching for New Physics. However, due mainly to the pollution of non- 
perturbative effects of strong interactions, an order of 1% CP asymmetries in D meson 
decay are still understandable in SM [19, 20, 21]. To date, no CP violation in charm sector 
is established. Though there were some hints of CP violation for the channels D ^ tttt and 
D KK [22, 23], the latest result from LHCb collaboration showed however, no evidence 
of CP violation in these channels [24]. 

A measurement of CP violation in the three-body decay has also been 

performed by LHCb [25]. With very high statistics, no localized or overall CP asymmetries 
are found. As has been shown in some aforementioned three-body decay channels of B 
meson, localized CP asymmetries in phase space can be enhanced by the interference of the 
decay amplitudes corresponding to two intermediate resonances with different spins, whereas 
the same mechanism should also apply to D meson decays. In this paper, we will perform 
a theoretical analysis of CP violation and differential decay width in the channel 
7r^7r+7r“, and will pay special attention to effects caused by the interference of intermediate 
state p^(770) and /o(500). 


2 Formalism of decay amplitudes and CP asymmetries 


In the region of the phase space around the resonances p^(770) and /o(500), the process 
TT^TT^TT is dominated by two cascade decays, 7r^/o(500) ^ TT^TT^TT and 

7r^p^(770) ^ 7r^7r+7r“. For the two weak decays 7r^/o(500) and 

7r^p^(770), the corresponding effective Hamiltonian can be expressed as [26, 27] 




KgVg(ci0? + C20|) 


i+h.c., 

i=3 ) 


( 1 ) 


where Gf is the Fermi constant, Vq-^q^ {qi and q 2 represent quarks) is the CKM matrix 
element, (z = 1, • • • , 6) is the Wilson coefficient, and Oi is the four quark operator, which 
can be written as 

01 = - 75)5/39/37''(1 - 75)Ca , 

O 2 = ^7a*( 1 - 75)557'‘(1 - 75)c , 

O3 = ujuil - 75)0'^ q'7^{I - 75 ) 9 ' , 

gf 

O 4 = «a7A*(l - 75)C/3 9^7''(1 - 75)9a , 

q' 

O5 = M7^(l - 75)cy^gV(l+ 75)9' , 

q' 

0 & = ^a7A*(l-75 )c/ 3 5^9^37'‘(1+ 75)9a , (2) 

q' 
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with Of and Of being tree operators, O 3 — Oq being QCD penguin operators, a and /3 being 
color indices, and running through all the light flavour quarks. 

The effective Hamiltonian for the strong decays tt+tt” and /o(500) ^ tt+tt” can be 
formally expressed as 


■^pOTTTr = - TT+d^^TT ), (3) 

^/OTTTT = 5/o7r7r/o(27r+7r“+7r°7r°), (4) 

where fo and tt^ are the field operators for p^, /o(500) and tt mesons, respectively, 
and p/oTTTT ^1*6 the effective coupling constants, which can be expressed in terms of the decay 
widths as 


487r 


"Jp'K'K 


fl’/oTTTT 


3/2 


^/o —)^7r+TT “ 


■ )-7r+TT 

rUn 
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Both /o(500) and p^(770) decay into one pion pair dominantly through strong interaction, 
and the isospin symmetry of the strong interaction tells us that T^o ~ VpQ ^^+^-, and T ^ 
— T 

2-*- /o^7r+7r“ • 

The decay amplitudes for the cascade decays 7r+/o(500) ^ 7v^7v^7v~ and 

7r+p^(770) ^ 7 T^ 7 T^ 7 T~ take the form 


( 5 low5high) 


•^/o low 5 high) — 


^/2GF^p7r7r^p('5high ^^low) 
'^low Tnp2 + iuipYp 


a4 — 


{VudV^d 


'V2 


aifpFi + a2fnAo 


2mlaQ 


FubFfjj ■ , 

LV {mc + md){mu + md) 
\/2GF{m% - ml)gfQ„„f^Fo 


firAo 
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■Slow ^/o F ^^/qT'/o 


-Vubv:, 
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2mlaG 
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• {KdV>2 

}. 
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respectively, where 5iow and Shigh are the invariant mass squared of tt+tt pairs with lower 
and higher invariant masses, respectively, = (5high,max + 5high,min)/2, with 5high^jnax(min) 

being the maximum (minimum) value of Shigh allowed by phase space for each siow, Fq^ Fi and 

Aq are short for the form factors and ^ respectively. 

All the a^’s are built up from the Wilson coefficients q’s, and take the form = q + Ci^i/Nc 
for odd i and = q + Ci-ijNc for even i. 

In the phase space that we are considering, the total decay amplitude for 7 t^ 7 t^ 7 t~ 

is dominated by M/q and M.pO. As a result, it can be expressed as 


M 


('^low5'^high) “1“ AdpO (siow?'^high)^ “1“ [-^low ^ -^high]? 


(9) 
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where 6 is the relative strong phase between the two amplitudes and Ad^o, which in 

principle, arises from long distance effect, [siow ^ ^high] represents a term which take the 
same form as that in the first square bracket except an interchange between siow and Shigh- 
For the calculation of the decay amplitude of the CP conjugate process D~ 7r+7r“7r“, 
which will be denoted as Ad, all one need to do is to replace the CKM matrix elements in 
Ad with their complex conjugates. 

The differential CP asymmetry for is defined as 


|MP + |Ai|2’ 

while the localized CP asymmetry can be expressed as 

d^highd^iowdAd I — |Ad| ) 
/^dshighC^siowd-MP + I-Mp)’ 

where, R represents certain region of the phase space that we are considering. 


( 10 ) 


( 11 ) 


3 Numerical analysis 


Table 1 list the input parameters and the corresponding references we used in this paper. In 
the following, we give some comments on these input parameters. We use the Wolfenstein 
parameterization for the CKM matrix elements, which up to the order of , can be expressed 
as [28, 29] 

Vud = 1 —2 - ^ ~ ~ ~ 

Vcd — “ 2(p + ir?)] + —^^(/) + i?7), 

Vub = \^A{p-iri), 

Vet = A\^ + rf), ( 12 ) 

with A, p, 77 , and A being the Wolfenstein parameters. To all orders in A, the relation between 
p, 7] and p, p can be expressed as [29] 


Vl — A‘^X^{p + if]) 

^ ^ Vl — A^[l — A^A'^(p + if])] 


(13) 


For the invariant mass dependence of the form factors p^^'^ and we use a model 

from Ref. [30], which take the form 


F{s) 


m 


777 , „ ^ 777 , „ ^ 


(14) 
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Table 1: Input parameters used in this paper. 


Parameters 

Input data 

References 

Fermi constant (in GeV“^) 

Gf = 1.16638 X 10"^ 

[32] 

Wilson coefficients 

Cl = -0.6941, C 2 = 1.3777, 

C 3 = 0.0652, C 4 = -0.0627, 

C 5 = 0.0206, ce — —0.1355, 

[27] 

Masses and decay widths 
(in GeV) 

mjj± = 1.86961, rjj± = 1.040 x 10“^^5 

BR{D^ 7T~^7T~7T~^) = 3.18 X 10“^ 

^p0(770) = 0.775, F^O(77o) = 0.150, 

^/o(500) = 0-5, ^fo{500) = 0-5, 

ttItp = 0.13957, 

rriu = 0.0023, rrid = 0.0048, 

rus = 0.095, rUc — 1.275, 

[32] 

Form factors 

= 0.67, Aq^P{0) = 0.64, 
a'^ = 1.19, = 0.36 

aP = 1.07, bP = 0.54 

[30] 

= 0.33, 

- 

Decay constants 
(in GeV) 

A = 0.13041, fK = 0.1562, 

[32] 

fp = 0.216, 

[33] 

Wolfenstein parameters 
of CKM matrix 

A = 0.22548+;;;[;K[;^t ^ = 0.810+K;[;^t 

P = 0.145W, 17 = 0.3431H11 

[34] 
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where F and X can be and tt, or and p, respectively. The form factor 

F^^f^{m^) which we use here is a rough estimation, and is consistent with branching ratio 
of /o(500)7r+ extracted from Dalitz analysis of Data [31]. 

As is observed by LHCb, the CP asymmetries around the vicinities of /o(500) and p^(770) 
have opposite signs for small and large values of Shigh in the case of B meson decay channel 
7 r^ 7 r+ 7 r“ [14]. In view of the above, for the case of we will focus on 

CP asymmetries of two regions, denoted as D+ and where D+ (D“) represents phase 
space satisfying Shigh > (<)^siow’ -^high > The CP asymmetry difference of the 

aforementioned two regions is 

AAcp = Aqp — A^p. (15) 

In Fig. 1 , the CP asymmetries A0p and their difference AAcp are shown as a function 
of the strong phase 5, where the strong phase 6 is assumed as a constant with respect to 
Slow nnd Shigh- It can be seen from Fig. 1 that AAcp is negative(positive) when 6 is around 
0 (tt). The magnitude of AAqp can reach as large as 0.5 x 10 “^ for some values of 8 . 
Especially, Fig. 1 shows that our mechanism indicates possibilities for AAcp being zero. 
This is interesting because that the experimental result from LHCb collaboration shows no 
CP violation in this channel. One can read off two zero points for AAqp from Fig. 1 , which 
are 8 i = 4.50 and 82 = 1.06. 

Figures 2 and 3 present in the phase space for 8 — 4.50 and 8 = 1.06, respectively, the 
relative differential decay width 7 of D+ ^ 7 r+ 7 r“ 7 r+, which is defined as 


T('^1ow 5 '^high) — 




Fp+_^^+^ — ^'^low^'^high 2567r^77lpF 


I-Mp 


(16) 


For comparison, we also present the relative differential decay width in Fig. 4 for 5 = 0, and 
that in Fig. 5 for the situation only resonance p^(770) is taken into account. Experimental 
data from LHCb shows that symmetries of event distribution around the p^(770) resonance 
are badly destroyed. The number of events around the p^(770) resonance for Shigh < ^siow 
are much larger than that for Shigh > shown in Ref. [25]. Besides, LHCb results 

also shows an enhancement of event distributions in the region of phase space where ^/s\^ 
and y/5high are around the masses of /o(500) and p^(770), respectively. These behaviours are 
roughly the same as those shown in Fig. 2, which indicate that our mechanism is consistent 
with experimental data when 8 = 4.50. 


4 Discussion 

We used a naive factorization approach for the weak decay processes in the calculation of the 
decay amplitudes. The reason is simply because that large part of the region in phase space 
that we focused on is off the mass shells of p^(770) and /o(500), the advantages of factorization 
approaches such as perturbative QCD [35], QCD factorization [36], Soft Collinear Effective 
Theory [37] are smeared out by the off shell effect. 

In determining the strong phase 5, we used the CP asymmetry difference of two regions 
of phase space instead of the differential CP asymmetry. The reason is because that the 
use of the differential CP asymmetry as a tool to determine the strong phase 8 is not an 
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appropriate approach at all. On one hand, if use the differential CP asymmetry, one would 
find that no strong phase can accommodate that with the data. On the other hand, if one 
check the nonzero differential CP asymmetries distributed in phase space for S = 4.50, one 
would find that the large differential CP asymmetries goes always with very small differential 
decay amplitude A4, indicating a cancellation between and A4 /q. In this situation, 

the dominance of these two amplitudes is no longer valid. Consequently, in order to deduce 
the differential CP asymmetries in this kind of regions, we should in principle consider other 
contributions to the decay amplitude A4 in these phase space regions, which is out of the 
scope of this paper. 

We choose the right boundary of the two regions and 14“ in phase space to be m^. 
Although it is not an unique choice, the boundary should not be far away from the vicinity 
of p^(770), in which case, the allowed strong phase 6 is not sensitive to the choice. On the 
other hand, either it is too large or too small than m^, the dominance of two resonance 
and /o(500) of the total decay amplitude is no longer valid. 

Other resonances can also contribute to the decay amplitude. For resonance such as 
/o(980), only a small part of the total resonance lies in the region. As a result, this con¬ 
tribution is small compared with p^(770). Resonance uu can enter in the amplitude through 
an isospin breaking effect, which is called — cj mixing mechanism. This mechanism can 
generate large differential CP asymmetries in the vicinity of 14. However, the width of uo is 
small, its contribution to regional CP violation is small. More importantly, the contribution 
oi u to CP asymmetry is independent of Shigh, and hence has no contribution to AAcp 

From Fig. 1, one can see that the CP asymmetries of the two regions 14+ and 14“ are 
nonzero for 6 = 4.50, which seems to be a disadvantage of our work. However, these two CP 
asymmetries are small, which both are 1.1 x 10“^. In principle, these small CP asymmetries 
are understandable by, for example, an inclusion of siow or Shigh dependence on 6. What 
important is, our division of the phase space /o(500) and p^(770) enlarge the effect of CP 
violations caused by the interference, consequently, can be used to determine the strong 
phase. 

5 Conclusion 

In this paper, we study the localized CP violation and differential decay width of the decay 
channel 44+ ^ 7r+7r“7r+. We focus our attention on phase space where the invariant mass 
of 7r+7r“ are around the vicinities of p^(770) and /o(500). We consider a mechanism which 
can generate large CP asymmetries on three-body decays of B meson, that is localized 
CP asymmetries caused by the interference of amplitudes corresponding to resonances with 
different spins. We found that by properly choosing a relative strong phase 5, the interference 
of amplitude corresponding to resonances /o(500) and p^(770) gives predictions that are 
consistent with experimental data both on CP asymmetries and differential decay widths. 

^Strictly speaking, since there are two identical particles in the final state for the terms 

of amplitudes are doubled by [siow ^ Shigh], as is shown in Eq. 9. As a result, the contribution of resonance 
uj in AAcp cannot be cancelled exactly. Besides, the interference between the amplitudes corresponding to 
/o(500) and lu is also small due to the smallness of Ct;’s width. 
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Our results generate no CP asymmetry differences {AAcp = 0) when the strong phase 
6 — 4.50. In the same time, the behaviour of event distribution around the vicinity of 
p^(770) and /o(500) is also understandable. 
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Figure 1: The CP asymmetry difference AAqp (solid line), the CP asymmetry of region 
(dashed line) and (dash-dotted line) as a function of the strong phase 6. 
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Figure 2: Differential branching ratio (in GeV of D~^ —>■ 7r"*‘7r tt"*" when (5 = 4.50. 
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Figure 3: Differential branching ratio (in GeV of D~^ —>■ 7r"*‘7r tt"*" when (5 = 1.06. 
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Figure 4: Differential branching ratio (in GeV of D~^ —>■ 7r"*‘7r tt"*" when (5 = 0. 


14 





(GeV) 



Figure 5: Differential branching ratio (in GeV of tt^tt tt^, where only the ampli¬ 

tudes corresponding to p^(770) are included. 
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